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ABSTRACT 

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range 
z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures 
yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time 
compared to previous HST imaging in this field, we probe ~ 0.65 (0.25) mag fainter in absolute UV 
magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the 
galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses 
or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that 
balances high redshift completeness and a low expected contamination fraction. We improve the 
number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter 
archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute 
magnitude limit of A/tjv = —17 and find a faint end Schechter slope of a = — 1.87lto 17- Using a 
similar color-color selection at z ~8 that takes account of our newly-added imaging in the F140W 
filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust 
estimate of the faint end slope at z ~8, a — — 1.94lg J,^. We briefly discuss our results in the context 
of earlier work and that derived using the same UDF12 data but with an independent photometric 
redshift technique (McLure et al 2012). 

Subject headings: cosmology: reionization — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Great progress has been made in recent years in stud- 
ies of the population of star-forming galaxies at redshifts 
z ~ 7 — 8. Following installation of the infrared Wide 
Field Camera 3 (WFC3), on the Hubble Space Tele- 
scope ( HST), the number of candidates has risen from 
a fe w dBouwens et alj I2008D to ~ 100 (IMcLure et al " 



20lllBouwens et aLH2010bl: IQesch et a' . 2010l|Yan et a' 



20101 ). In addition to providing hints of the early galaxy 
population to z ~ 8, previous data sensitive to z ~ 7 
galaxies have provided initial determina tions of their 
rest-frame UV colors, stellar populations (IMcLure et al.1 
l20lHlBouwens et al.ll2012bHDunlop et al ll2012af). stellar 
mass e s and likely ag es (Labb eet alj2010tl Gonzalc z et al] 



120101 : iMcLure et all 120111 : iFinkelstein et al 



nebular emission line strengths (jLabbe et al.l 



2012) 



2012) 



and 
Our 

work builds upon these previous efforts to present the 
first drop-out selected samples and luminosity function 
determinations for redshift z ~ 7 and z ~ 8 sources 
from the 2012 Hubble Ultra Deep Field project (here- 
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after UDF12; GO 12498, PI: R. Ellis). 

Before UDF12, progress has naturally been greatest 
at redshift z ~ 7 where synergy between ground- and 
space-based surveys has effectively exploited the full dy- 
namic range of acc essible galaxy lum inosities. Early sur- 
veys from Subaru (lOuchi et aljl2009h. and the ESO Very 
Large Telescope (jCastellano et nlJ"[2T)l()!. have probed 
the luminous component of the star-forming population 
over an area > 1000 arcmin 2 . More recently, the Ul- 
traVISTA survey has covered 3600 arcmin 2 in the COS- 
MOS field, locating z ~ 7 galaxies to M uv = -22.7 
(Bo wler et al.ll2012f ). 

However, only HST can probe the important faint end 
of the galaxy luminosity function at these redshifts. An 
early result from the 2009 Hubble Ultra Deep Field cam- 
paign (GO 11563, PI: Illingworth, hereafter UDF09) was 
the discovery of a n abundant popula t ion of sub- luminous 
galaxies at z c± 7 (IQesch et alJl20lfl lBunker et"aTl l201Ct 
IMcLure et al.ll2010f ) corresponding to a Schechter faint 
end slope a between -1.7 to -2.0. In such a distribution, 
the bulk of the integrated luminosity density arises from 
low luminosity galaxies that may be responsible to m ain- 
taining cosmic reionization (Ro bertson et al.ll2010aft . 

Clearly the luminosity function of star-forming galax- 
ies at redshifts z ~ 7-8 is of great importance. How- 
ever, given the large range in luminosity that must be 
sampled, wider-field HST surveys have proved an im- 
portant complement to panoramic ground-based surveys. 
WFC3 data from the GOOD S Early Release Science 
ERS1 rt Windhorst et al.l [20T1. and CAND ELS fields 
Grogin et al.l 1201 It iKoekemoer et al.l 120 111 ) have sam- 
pled intermediate luminosities — 21 < Mtjv < — 19. The 
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HIPPIES and BoRG pure parallel su rveys have pro- 
vided additional candidates at z ~ 8 (|Yan et al.l 120111 : 
iTrenti et al.l[20TTI: iBradiev et al.ll2012T ). 

Collectively these surveys have provided a reasonably 
clear view of the galaxy luminosity function at z ~7 at 
the luminous end but there remain disagreements regard- 
ing the precision with which the faint end slope is deter- 
mined. While the UDF09 collaboration has measured a 
faint- end slope of a z ^ = —2.01 ± 0.21 (|Bouwens et al.l 
1 2 llf ) incorporating the UDF, parallel fields, and the 
ERS data, a competing determination utilizing the size- 
luminosity relation measured from the UDF09 data 
and the CANDELS Deep+Wide surveys in three fields 
finds a shallower fai nt-end slope of a z ^ = —1.7 ± 0.1 
(|Grazian et al.ll2012l ). The luminosity function at z ~8 is 
even more uncertain, both because of the limited depth of 
the necessary photometry (jDunlop et al.ll2012a[ ) and the 
possibility of contamination from lower redshift sources. 

Distant star-forming galaxies are commonly found us- 
ing sojn£VB£iant of the Lyman break technique pioneered 
bv ISteidel et al.l (|1996l ). At z > 6.5, the opacity due to 
neutral hydrogen in the intergalactic medium means the 
source flux below a rest- wavelength A r es t = 1216A is 
dimmed by factors > 5 (|Madaul I1995D . A search ex- 
ploiting this effect has utilized one of two methods. In 
the 'dropout' technique, objects are selected within a 
carefully-chosen window in color-color space specifically 
tuned to select star-forming galaxies within the required 
redshift range while minimizing the contribution from 
lower redshift contaminants. At z ~ 7 (where the Ly- 
man break falls near the overlap of the HST zgso and 
Y105 filters) the Lyman dropout is chosen via a red color 
in z — Y , and the star-forming nature of the galaxy via 
a blue color in Y — J. Additionally candidates are re- 
quired to lie below a certain threshold in deep optical 
data; this further limits contamination by lower redshift 
sources. Early demonstrations of th e drop-out technique 
at r edshifts z > 6 w e re pres ented by lBunker et al.l (2004) 
and lBouwens et al.l (|2004al ). 

An alternative approach uses the full array of 
broadband detections and upper limits in the con- 
text of photomet ric re dshift codes (IMcLure et al.l 120101 : 
iFinkelstein et al.ll2010t: iMcLure et al.ll201lD . These codes 
employ a range of synt hetic spectra for galaxies (e.g., 
IBruzual fc Charlotll200l over all redshifts of interest and 
optimum fits are produced for each source in the catalog. 
As in the dropout selection technique, leverage comes 
primarily from the Lyman break but the method is par- 
ticularly advantageous when detections are available in 
more than 3 filters. 

In general, a greement between th e two techniques is of- 
ten quite good (jMcLure et al.ll2011l ). However, the meth- 
ods are distinct and each has features relevant for inter- 
preting their photometric samples. The SED method 
provides redshift probability distributions for individual 
sources but may be susceptible to systematic errors in- 
herent in population synthesis models, such as uncer- 
tainties in the reddening law and star-formation histo- 
ries, when differentiating between low-redshift interlop- 
ers and true high-redshift sources. In contrast, our drop- 
out selection only requires the observed color informa- 
tion independent of stellar population synthesis model- 
ing, but requires careful simulations to quantify the pos- 



sible presence of low-redshift contaminants satisfying the 
break criterion. Given their complementarity and differ- 
ing systematics, independent luminosity function deter- 
minations from both methods will be helpful in further- 
ing progress. The goal of the present paper is there- 
fore to exploit this unique data to provide the best cur- 
rent constraints on the UV luminosity function of star 
forming galaxies at redshift 7 and 8 usin g the dropout 
techn ique. A companion UDF12 paper ([McLure et al.l 
I2012T ) presents the results of a search through the data 
for z > 6.5 candidates using the photometric redshift 
technique. 

The present paper is one in a series devoted to scien- 
tific results from the UDF12 campaign, which provides a 
significant advance over the earlier UDF09 imaging in 
this field appropriate for the present luminosity func- 
tion study and the c omplementary analysis discussed in 
IMcLure et all (|2012t ). The UD F 12 survey design and it s 
data processing is discussed in iKoekemoer et al.l (|2012|) . 
Public versions of the final reduced WFC3/IR UDF 12 
images, incorporating all earlier UDF data, are avail- 
able to the community on the UDF12 web pag^3- Initial 
z > 8.5 detection s in the UDF 12 data were presented by 
lEllis et all ()2012f ). while th e UV continua of high- redshift 
sources were measured by iDunlop et al.l (|2012bl ) . A re- 
view of the implications of the survey in the context of 
cosmic reionization including the results of the present 
paper is provided in iRobertson et al.1 (|2012t ) . 

A plan of the paper follows. In Section 2 we introduce 
the UDF12 data and the brighter ground-based and HST 
data sets and their reductions essential for realizing an 
analysis of the luminosity function at z ~7 and 8. Sec- 
tion 3 discusses the important decisions we have taken in 
color-color space to optimize the completeness of galaxies 
at these redshifts, while minimizing contamination from 
lower redshift sources. Section 3 also presents the final 
list of galaxies in the two redshift intervals that we use 
for our analysis. In Section 4 we present our luminosity 
functions and in Section 5 we briefly discuss our results 
in the context of earlier work, highlighting the important 
advances made possible through the UDF12 campaign. 

Throughout this paper, we adopt a A— dominated, 
flat universe with 17a = 0.7, flju — 0.3 and Ho = 
70h.7o kms -1 Mpc -1 . All magnitudes in this paper are 
quoted in the AB system (|Oke &: Gunnl 119831 ). We will 
refer to the HST ACS and WFC3/IR filters F435W, 
F600LP, F606W, F775W, F814W, F850LP, F105W, 
F098M, F125W, F140W, and F160W as B 435 , V 600 , V 606 , 
«775, «8i4, ^850, ^098, ^105, ^125, ^wo, and H w0 , respec- 
tively. 

2. DATA 

Central to any analysis of the galaxy luminosity func- 
tion is the collation of a complete sample of galaxies 
within the chosen redshift interval spanning a wide range 
in luminosities, free from bias and with any interlopers 
minimized. In this section we introduce the UDF12 and 
more luminous auxiliary datasets. 

2.1. UDF 

To provide the best constraints on the faint end of the 
luminosity function at z ~ 7 and 8, the primary ad- 

8 |http://udfl2 .arizona. edu/ 
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vance of this paper is the increased depth and redshift 
fidelity provided by our new UDF 12 su rvey. The UDF 12 
progr am dataset (described in full in iKoekemoer et all 
120121 ) represents a significant improvement over the pre- 
vious UDF09 observations in several respects. In par- 
ticular, the survey was purposefully designed to improve 
our understanding of the redshift z ~ 7 and 8 luminosity 
functions. Firstly, we increased the total exposure time 
in the key Y105 filter over that in the UDF09 campaign 
by a factor of 4, with the addition of 71 new orbits. As 
this filter lies at the edge of the Lyman break for galaxies 
at z ~ 7, this ensures we can probe significantly fainter 
(by ~0.4 mag) in absolute magnitude at z ~7 and with 
considerably grea ter fidelity in redshift selection at z ~8 
([Ellis et al.ll2012f > . A further improvement is the addition 
of comparable imaging in the newly-utilized J140 filter. 
By stacking our detections in this filter with either those 
at J125 at z ~7 or H\qq at z ~8, we secure improved 
detections that correspond to extending the depth by an 
additional ~0.1 mag in each case. 

Our fina l analysis is based on the compilation dis- 
cussed by Koeke moer et all (|2012D which incorporates 
all earlier WFC3 imaging in the UDF including the ear- 
lier UDF09 campaign (GO 11563; PI: IUingworth) and 
less deep imaging undertaken as part of the CANDELS 
survey (GO 12060, Pis: Faber & Ferguson). In to- 
tal the imaging constitutes 100, 39, 30, and 84 orbits 
in each o f the Y lnR , ■7i2K. Ji4n a nd Hi 60 filters respec- 
tively (see lKoekemoer et al.ll2012l for further details). An 
important associated dataset in this field is the ultra 
deep ACS imaging da ta from the 2004 UDF campaign 
(|Beckwith et al.ir2006f) essential for a further rejection of 
low redshift sources. 

2.2. Auxiliary data 

To constrain the bright end of the z ~ 7 and 8 lumi- 
nosity functions we take advantage of several auxiliary 
WFC3 datasets which are somewhat shallower than our 
UDF12 data but nonetheless unique in their coverage 
and depth. At z ~ 7, we include the UDF parallel fields, 
UDF-P1 and P2 (sometimes referred to as UDF-P12 
and P34, respectively) from various surveys including 
UDF05 (GO 10532, PI: Stiavelli) and the aforementioned 
UDF09 survey, the Early Release Science(ERS) WFC3 
campaign from the WFC3 Science T eam (GO 11359, 
PI: O'Connell: [Windhorst et 1012011 . and imag ing m 
the GOODS-S field from CAN DELS-Deep (jGrogin et al.l 
120111: IKoekemoer et all 120111 ). We also ad opt the dat- 
apoint s from the ground-base d su rveys of lOuchi et all 
(20091 . [Castellano et all pOlOh and lBowler et al.l (|2012j ), 
which are provide vital observations of the rare popula- 
tion of galaxies brighter than Muv,*- 

At z ~ 8, it is difficult for ground-based pro- 
grams to provide any meaningful constraints, so we in- 
stead include data provided by two HST pure-parallel 
WFC3 program s: HIPPIES (GO/PAR 11702, PI: Yan; 
Yan et al.ll20Tl and BoRG (GO /PAR 11700, PI; Trenti 



Trenti et all 1201 It iBradlev et ail 120 121 ). As the shortest 



fully identify galaxies at z ~ 8, and therefore we do not 
use them for our selection at z ~ 7. 

We summarize the filter coverage, survey area, limiting 
depths in the selection filter(s) and the number of high 
redshift candidates for each of these auxiliary datasets in 
Table 1. 

2.3. Data Reduction 

Prior to applying photometric color cuts optimized for 
the selection of z ~7 and 8 galaxies, each survey data 
was similarly reduced to provide a series of processed 
and calibrated WFC3 frames. We describe below the 
data reduction steps taken for each field. 

2.3.1. UDF and Parallels 

The preliminary processing stages that yield images 
ready for source select ion are discussed in detail by 
IKoekemoer et all (|2012t ). Briefly, we first process the 
raw images using the Pyraf/STDAS task calwf3, which 
flags bad pixels, and corrects for bias and dark current 
throughout the detector. After processing, the images 
are then registered and stacked using a versio n of the 
MultiDrizzle algorithm (Koeke moer et al.l 120111 ) to cre- 
ate the final mosaics on a pixel scale of 0.03" per pixel. 
This processing was carried out to create reductions of 
the UDF, UDF-P1, and UDF-P2 fields. 

2.3.2. ERS and CANDELS-Deep 

For the ACS data in the CANDELS-Deep and ERS 
fields, we use the pu blicly available v2.0 m osaics from the 
GOODS campaign ( Gi avalisco et alll2004D . We augment 
this with our own compilations of both the isi4 and Zsso 
data taken in parallel during the CANDELS WFC3 cam- 
paign. We combine the publicly available single epoch 
mosaics for these filters weighting by expos ure time, us- 
ing th e image combination routine SWARP (|Bertin et all 
2002). In the case of the zsso data, we add this to the 
already existing GOODS mosaic. 

We combi ne the public WFC3 mo saics from the CAN- 
DELS team (jKoekemoer et alll20lil ) in the same manner 
across CANDELS-Deep. For the ERS WFC3 da ta, we 
use the reduction described in lMcLure et all (|201 lh . Due 
to the lack of sub-pixel dithering available in the wider 
area fields, the CANDELS and ERS mosaics were pro- 
duced with final pixel scales of 0.06" per pixel. 

2.3.3. BoRG + HIPPIES 

Data from the BoRG and HIPPIES programs differ 
significantly from the prior data sets, as these are pure 
parallel HST sur veys. Details o f the obser vation strategy 
can b e found in iTrenti et all (|2011h and I Bradlev et al l 
2012). The complete dataset presented in IBradlev et all 
20121 ) consists of 59 independent fields. Due to the dif- 



wavelength coverage of these surveys is provided only 
with the Vqoqw (or V eQ0L p), Y 098 , and J125 filters, a ro- 
bust spectral break can only be verified between Yogs and 
J125, as the wavelength spanned between V6O6W and Yogs 
is too great and there remains no optical rejection filter 
at shorter wavelengths. Thus, these fields can only use- 



ferent nature of dealing with pure parallel data of lim- 
ited depth, and the survey's high expected contamina- 
tion rate of ~ 40%, we make no effort to undertake an 
a nalysis ourselves. In stead, we adopt the data points 
of IBradlev et al.l (|2012l ). and in Section 4, determine fits 
with and without these data added. 

2.4. Photometry 

In the case of the UDF12, which pushes HST to new 
limiting depths, we adopted a robust technique to locate 
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and measure the fluxes from each faint source, rather 
than relying on e rrors output from the S Extractor source 
extraction code (jCasertano et al.l [2000f ). As the corre- 
lated noise produced as a result of applying the Mul- 
tidrizzle algorithm produces a subtle underestimate of 
uncertainties when using SExtractor, in the specific case 
of the UDF12 we chose to use our own IDL photometry 
routine to compute all fluxes quoted in this paper. 

We briefly summarize the various stages. Processing 
proceeds by using a % 2 stack of all the images to iden- 
tify regions of blank sky over the area in question. A 
grid of blank apertures is then generated, with separa- 
tion larger than the pixf rac footprint of Multidrizzle, 
ensuring that noise between adjacent apertures is not 
correlated. To estimate the level of any residual back- 
ground around an object of interest, we take the median 
flux of the 50 closest blank apertures. We adopt the 
root mean square of the flux in these blank apertures as 
our uncertainty in the background level. Flux es are then 
comp uted using the APER routine in IDL ([Landsman! 
1993). For the shallower non-UDF12 fields, although we 
utilized SExtractor to compute fluxes and background 
levels, the flux uncertainties were still estimated using 
this improved technique. 

As the HST point spread function (PSF) varies signif- 
icantly with wavelength (particularly between the ACS 
optical images and near-infrared images from WFC3), it 
is important to account for this change when measuring 
accurate colors. To improve detections and color mea- 
surements for the faintest sources, aperture sizes should 
be quite small yet properly account for wavelength- 
dependent PSF variations. We chose circular apertures 
whose diameters encircle > 70% of the total flux from 
a point source. For the ACS -B435, Vqqq, 1775, isu, and 
Z850 filters, we adopt aperture diameters of 0.30" in all 
fields processed with 0.03" pixel diameter (UDF12, UDF- 
Pl, and UDF-P2), and 0.40" in all other fields. For the 
WFC3 filters, the aperture diameters are 0.40", 0.44", 
0.47", and 0.50", for Y 105 /Y 098 , J 125 , J U o, and H 160 , re- 
spectively. Such small apertures remain meaningful both 
because of the FWHM of the HST PSF (ranging from 
0.09" in B435 to 0.17" in H 160 ), and the precise degree 
of alignment of the individu al image sub-exposures (bet- 
ter the 0.005" i n the UDF ; jKoekemoer et al.l 120121). In 
a related paper, IQno et al.l (|2012f7 validate these choices 
by measuring half-light radii, r/,.;, for stacks of our high- 
redshift samples, and find values of 0.36 kpc at both 
z ~ 7 and 8, which translate to angular half-light radii 
of 0.07" (not including any broadening due to the PSF). 

3. CANDIDATE SELECTION 

We now turn to the photometric selection of star- 
forming galaxies at z ~ 7 and 8, using the photometric 
catalogs generated as discussed in Section 2. A key is- 
sue for our dropout selection technique is the optimum 
choice of the two color cuts used to select candidates. 
The goal is to balance completeness in high z selection 
against a low fraction of foreground interlopers and spu- 
rious sources. As the available filter sets differ for each 
of the component surveys, we discuss each case in turn. 

3.1. Potential contaminants 

A selection criterion tuned to select only high-redshift 
galaxies at high confidence, with essentially no contami- 



nation, would be impractical leading to a severely limited 
sample. In order to be inclusive in such a search, a crucial 
condition for an accurate determination of the luminosity 
function, a small degree of contamination can be toler- 
ated provided the fraction is reasonably well-understood. 

In fact, the primary sources of contamination in high- 
redshift searches are now well-known. We briefly review 
them here. Low-temperature galactic dwarf stars dis- 
play quite similar c olors to high redsh ift galaxies in the 
near infrared (e.g.. iBowler et aT1l2012| ). While ground- 
based surveys can only constrain stellar contamination 
by comparing SED s of cool stars to observed colors 
(jBowler et al.ll2012l ). WFC3 data has the distinct advan- 
tage of an extremely sharp point-spread function, rang- 
ing from 0.15" in Y105 to 0.18" in H160. Previous studies 
in the UDF have utilized this to conclude that Galactic 
stars are not expected to contribute significantly ( < 2%) 
to contamination in extremely deep, small a rea surveys 
(jFinkelstein et al.ll2010HBouwens et al.ll201lD . 

A more relevant concern is potential interlopers at 
z ~ 2. Around this redshift, the Balmer and 4000A 
breaks will lie near the same observed wavelengths where 
we search for a Lyman break in our high-redshift sources. 
However, unlike the Lyman break, the Balmer break has 
a maximum possible depth, aiding us in isolating rob ust 
high-redshift sources (see Figure 2 of lKriek et alllMol for 
actual measurements at z > 0.5). Due to our extraordi- 
narily deep optical data, these objects must either then 
be severely reddened in the rest-frame ultraviolet or sig- 
nificantly affected by photometric scatter in order to be 
picked up by our selection criteria. We account for the 
possibility of such sources scattering into our selection 
window with our contamination simulations described 
below. 

3.2. Optical non-detection criteria 

Applying a rigorous set of optical non-detection crite- 
ria is key to obtaining a clean sample of high-redshift 
sources by removing the lower redshift contaminants we 
describe above, as we expect all intrinsic flux from true 
high-redshift sources to be nearly extinguished by neutral 
Hydrogen at these wavelengths (Madau 1995). Here we 
adopt a slightly mod ified version of the criteria used in 
iBouwens et al.l([201lD to eliminate sources with marginal 
optical detections from our selections. 

The criteria we apply are as follows: (1) The measured 
flux in each filter shortward of the dropout filter is less 
than 2.0cr. For the z-drops, this includes -B435, Vqqq, and 
1775; for the F-drops, we add zsso- (2) No more than 
one of the filters listed above shows a detection above 
1.5cr. (3) To effectively add all the optical data, we fi- 
nally require Xop*@ must be less than a threshold value, 
depending on the observed magnitude of the source. We 
compute this value using both our standard 0.30" diam- 
eter circular apertures, and a smaller 0.21" aperture, to 
rule out the most compact contaminants. At or below 
the 5cr limit of tiiab = 29.5 (in uncorrected, aperture 
magnitudes), we adopt a \opt upper limit of 2.5, while at 
the 10cr limit of tuab = 28.7, we relax this to limit 5.0. 

9 defined as Xopt = ^- , i s ^9 n (fi)(fi/ rJ i) 2 > where the i index runs 
across -B435, V606, and 1775 for z-drops; additionally «814 where 
available, and Z850 for y-drops 
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A linear interpolation is used to determine the limit for 
magnitudes between the 5 and 10cr level. 

3.3. Contamination simulations and the adopted 
UDF12 color-color selection 

To estimate the number of contaminants we expect, 
we utilized the excellent HST photometry for objects at 
25.0 < l/i6o < 27.0 (as measured in our 0.5" diameter 
apertures) in the UDF. As we can robustly rule out the 
presence of a Balmer break at this depth, we selected, as 
our base color distribution of contaminants, all sources in 
this magnitude range that displayed at least one optical 
detection. The key assumption in our simulations is that 
the color distribution of these potential contaminants is 
unchanged as one moves down to the fainter magnitudes, 
where the majority of our dropout galaxies lie. Consider- 
ing the trend in color in this magnitude range, we believe 
this may even be a conservative assumption. For our 
sample at 25.00 < Hi 60 < 25.25, the median z S50 — Yi05> 
Y105 — J\25i and J125 — i?i60 colors are 0.45, 0.31, and 
0.31, respectively. At 26.75 < H wo < 27.00, these three 
colors show even less extreme median colors of 0.32, 0.17, 
and 0.15. 

We then create an array of synthetic sources, matched 
to the actual number of observed sources in the UDF 
in bins of 0.1 magnitudes. To get an accurate rep- 
resentation of sources intrinsically below our detection 
limit that have some chance of being scattered upwards 
into detection, we extrapolate the number counts beyond 
iJi6o = 28.7 (equivalent to a source at 10. Oct) as a func- 
tion of magnitude via a power law down to sources in- 
trinsically as faint as l.Oer. The colors of these synthetic 
sources are chosen to obey the same color distribution 
as the brighter contaminant population described above, 
but noise consistent with their synthetic magnitudes is 
then added. These sources are then subject to the same 
optical non-detection criteria described above, but the 
colors cuts in z 850 - Y W5 and Yios - J125 (Yi.05 - J125 and 
J125 — Hieo) for z-drops (F-drops) are varied in steps of 
0.05 from as shown in Figure [T] To create these plots, 
we repeated the simulations for N = 1000 UDF fields. 

Clearly, the most robust constraint on eliminating con- 
taminants comes from requiring a large break in the bluer 
color for each sample, but gains are also made by lim- 
iting the color longward of the break to relatively blue 
values. Provided with these estimates, we select chose to 
implement the color selection criteria as shown in Figure 
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This leads to a selection for z-drops broadly within 
the redshift range 6.2 < z < 7.3. For the z ~8 study, 
we use our new ultra-deep WFC3 F105 as the dropout 
filter which leads to a sample spanning the redshift range 



7.3 < z < 8.5. A discu ssion of more dis tant sources in 
the UDF is provided bv lEllis et all ([2012I ). 

Finally, to ensure our sources are robust, we demand 
a detection significance of 3.5c in the filter immediately 
longward of the Lyman break (F105 for our z-drops, J125 
for our F-drops, and a similarly robust detection in one 
further WFC3 filter at longer wavelengths. 

These selection criteria are designed to provide as large 
a sample of galaxies as possible above redshift 6.5, while 
minimizing the effect of contamination. We plot selec- 
tion functions for both z-drops and F-drops in Figure 
2 (see Section 4.1 for details on the selection function 
simulation). At bright magnitudes J140 < 27.5, our z- 
drop color cuts provide a nearly complete census of star- 
forming galaxies between 6.30 < z < 7.15, while the 
F-drop cuts do the same between 7.35 < z < 8.60. How- 
ever, due to our strict optical non-detection criteria, and 
the area subtended by other sources in UDF12, our max- 
imum selection efficiency does not exceed ~ 65%. 

Extending our Y105 — J125 color cut for z-drops to 0.5 
would only add an extra ~ 0.1 in redshift space to our 
selection function, as the color tracks are rapidly depart- 
ing from our selection window as can be seen in Figure 
3. This would additionally add an extra 0.3 expected 
contaminants. Combined with the concern that our new 
F105 data is actually deeper than the existing z 850 data, 
and that the primary contaminants themselves are in- 
trinsically red, we opted for this conservative Y105 — J125 
color cut for our z-drop sample. Similarly adding an 
extra 0.1 mag to our J125 — i?i6o limit for F-drops is ex- 
pected to add ~ 0.4 contaminants. As can be seen from 
the F-drop color-color plot in Figure 3, the density of 
F-drops with J125 — i?i6o colors this red is quite low, so 
we chose to truncate the selection at J125 — -Hi60 < 0.4, 
to limit the contamination. Adopting these cuts and all 
of our previously discussed selection criteria, we expect 
2.79 and 1.42 contaminants per UDF field for z-drops 
and F-drops, respectively. 

We include candidate lists and photometry for our fi- 
nal selection of z-drops and F-drops in Tables 2 and 3, 
respectively, and color-color plots in Figure 3. In total, 
we select 47 z-drops, and 27 F-drops, of which 20 and 9 
were only identified using our new UDF12 data. 

3.4. UDF-P1 and UDF-P2 

The two UDF parallel fields, observed as part of the 
UDF05 (GO10632; PI: Stiavelli) and UDF09 (GO 11563; 
PI: Illingworth) campaigns comprise the two data sets 
most similar to our UDF 12 dataset, though the depths 
are ~0.5 mag shallower. As such, we utilize the same 
color-color criteria determined for our UDF12 selections. 
Because the optical data is shallower, we tighten our Xopt 
upper limit for selection to 1.5 (3.0) at the 5a (10a) 
aperture magnitude limit in each field. We utilize the 
ultradeep 128 orbit z 8 i4 ACS data taken in parallel to 
our main UDF12 program, which covers roughly 70 % 
of the P2 WFC3/IR field. For z-drops within this area 
we impose an additional criteria of z 8 i4 — F.05 > 2.0 OR 
/814/0814 < 2.0. 

3.5. EES 

For the ERS dataset, the F)98 filter was utilized, so 
we take care to alter our selection criteria accordingly. 
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We ch ose to use the criteria derived by iBouwens et al.l 
(|2011f ). at z ~ 7, which, despite the different filter selec- 
tion, produce a selection function that probes a similar 
range in redshift. These criteria are (1) zsso — ^125 > 
0.9, (2) z 850 - J 125 > 0.8 + l.l(Ji 25 - Fieo), and (3) 
zss n - J125 > 0-4 + l.l(F nfl 8 - Ji2h)- At z - 8, we adopt 
the lBouwens et all (|2011f ) Yogg — J125 lower limit of 1.25, 
but chose the same J125 — i?i60 < 0.4 cut we use for 
the UDF fields, to ensure the selection functions for our 
analysis are as homogeneous as possible. 

3.6. CANDELS 

For the CANDELS field, we use the same criteria as 
adopted for UDF-P34, including the addition of the zsi4 
criteria. Due to the varying depth of the data, we di- 
vide the CANDELS field into 3 distinct sub-regions when 
performing our analysis. The first region consists of 
the immediate area around the HUDF that is covered 
by ACS optical imag ing taken as part of the HUDF04 
(jBeckwith et alJfeOOGl ) campaign, but falls outside the re- 
gion probed by our UDF12 campaign, due to the smaller 
field of view of the WFC3 instrument. Here, the optical 
data is nearly a full 2 magnitudes deeper than the avail- 
able IR data, resulting in negligible contamination rates 
in our sample. 

The remaining sample is further divided into an East 
and West region. The East and West regions have iden- 
tical depth in the ACS data from the GOODS program, 
and in the J125 and H 160 filters, but the F105 depth in 
the West region is approximately 0.4 magnitudes shal- 
lower than that of the East. For the purposes of contam- 
ination and completeness simulations, we separate these 
two fields in order to properly account for the variation 
in depth. 

4. THE LUMINOSITY FUNCTION AT Z ~ 7 AND - 
8 FROM UDF 12 DATA 

With our candidate selection completed, we now turn 
to evaluating the z ~ 7 and 8 luminosity functions. The 
key issue in converting numbers of sources of known ab- 
solute magnitude into a comoving density of luminosi- 
ties is, of course, the redshift-dependent selection func- 
tion which defines the visibility as a function of apparent 
magnitude as well as the optimum algorithm for fitting 
a function to the resulting number density. We now in- 
troduce the algorithms we will utilize for both of these 
critical steps. 

4.1. Simulations 

We first describe how we calculate the selection func- 
tion used to determine the effective volumes for our lu- 
minosity function calculations. 

To create synthetic fluxes for galaxies in our sim- 
ulations, we assume an input UV slop e 8 = —2.0 
This ch oice is motivated by the results of iDunlop et al.l 
(2012b), who find no conclusive evidence for an intrinsic 
scatter in 8 from this value at z ~ 7 — 8. We compute 
fluxes applying intergalactic extinc tion from iMeiksin 
(2006) to a Bruzual and Chariot (|Bruzual fc C hariot 
2003) synthetic spectrum consistent with this value of 

We parameterize our selection efficiency, S(m,z), as 
a function of redshift and the magnitude for the filter 



(or filters for the UDF12 data) used to determine the 
rest-frame UV magnitude. In each field, we determine 
the shape of the selection function first using numeri- 
cal simulations only, which take into account the limit- 
ing magnitudes in each filter for point sources. In steps 
of 0.01 in redshift, and 0.05 in magnitude, we take the 
synthetic flux from our galaxy model, perturb it by the 
appropriate noise, and apply our color-color selection cri- 
teria. At each step of z, m, and for each field, we perform 
this N = 1000 times to construct a complete surface for 
our selection function. We define the selection function 
produced by this process as S numer i C (m, z). 

However, this selection function is only appropriate if 
these galaxies are point sources in otherwise blank fields, 
which is certainly not the case. The marginally resolved 
nature of our targets will result in higher incomplete- 
ness levels at faint magnitudes than for point sources. 
It is imperative to correct for this effect, as a varying 
completeness correcti on can produce larg e differences in 
the faint end slope (|Grazian et al.l [20121) . To account 
for this incompleteness, as well as that caused by area 
in the images obscured by brighter sources, we rely on 
inserting synthetic galaxy images into our mosaics for 
each field. We generate a synthetic template with a 
Sersic i ndex of 1.5, consis tent with a stack of LBGs at 
z ~ 4 (|Oesch et al.ll2010T) . The template image has a 
half-light radius of 0.35 kpc, motivated by the results of 
lOno et ahl (|2012l ) who perform a detailed measurement of 
the sizes of UDF12 high-redshift candidates. This tem- 
plate is then convolved with the point spread function 
for each filter, multiplied by the appropriate model flux 
as described above, and inserted into the image in a ran- 
dom location. After inserting N = 1000 non-overlapping 
sources, we run the SExtractor program for object detec- 
tion, and compute fluxes and errors in the same manner 
as for our science images. We perform this simulation 
at the peak of each of our selection functions in red- 
shift space, in steps of 0.05 in magnitude. This peak 
efficiency at a given magnitude e(m) is then used to nor- 
malize our selection function such that our total selection 
function used to compute effective volumes is given by 
Stot a i{m, z) = S numeric (m, z) x e(m). The final selection 
functions for both z-zdrops and Y-drops in the UDF are 
presented in Figure 2. 

4.2. Maximum Likelihood Luminosity Functions 

Using the new UDF12 data and previous observations, 
we assemble dropout samples at z ~ 7 and ~ 8 in mul- 
tiple fields. For each sample, we split the galaxy num- 
ber counts into bins of width AM = 0.5 in absolute 
magnitude Mf/y. We use these binned galaxy number 
counts and our simulations of photometric scattering and 
low-redshift contaminants to determine the high-redshift 
stepw ise maximum likelihoo d (SWML; Efstathio u et al.l 
1988) and Schcchter (1976) galaxy luminosity functions 
$(Muv) in units of Mpc -3 mag -1 . 

4.2.1. Stepwise Maximum Likelihood Luminosity Function 

The SWML luminosity function aims to jointly fit the 
binned galaxy abundance in the k th of iVbin magni- 
tude bins. For each bin, the maximum likelihood values 
for <S>k are determined by using the observed number of 
galaxies n b s ,fe, the effective volume V c g.k for galaxies 
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with intrinsic Myy in the bin, and the probabilities Pjj 
for scattering galaxies with intrinsic Muv in bin i into 
bin j owing to photometric noise, and the number of 
low-redshift contaminants n CO n,fe for the bin calculated 
as described in § section:simulations. Given we con- 
struct the expected number of galaxies in each bin as 



''cxp,/;; " 



k y c s,k 




+n con ,fe, (5) 

where the summations run over Nhm- In practice, when 
using a bin width AM UV = 0.5 galaxies do not scatter by 
more than one bin and the summations are trivial. We 
account for photometric scattering of sources into our 
faintest bin by a simple extrapolation of the luminosity 
function to yet fainter magnitudes. 

To fit the shape of our SWML, we can use the like- 
lihood of observing n b s ,fc given the expected number 

^-exp.fc • 



p(«ob s ,fc|n-cxp,fc) 



^cxp,/c 



(6) 



The SWML luminosity function is determined 
by maximizing the product of these individual likeli- 
hoods across all bins, and across all fields. Photomet- 
ric scatter correlates the individual $fe for each field, 
and the SWML values must therefore be determined 
simultaneously. This calculation amounts to an iVbi n - 
parameter estimation problem, and we use the MultiNest 
nested sampling code f or Bayesian infer ence problems 
(|Feroz fc Hobsonl I200& iFeroz et all [2001 to maximize 
the <&fc likelihoods. 

To fit the overall normalization of our SWML, we sum 
fiobs.fc across all fields for each magnitude bin. Since we 
expect this quantity to be Poisson distributed, we can 
easily generate a posterior distribution of the normaliza- 
tion for each bin. To find our final posterior distribution, 
for the normalization, we multiply the posteriors gener- 
ated in this manner across all bins. 

The results of the SWML luminosity function calcu- 
lation for redshift z ~ 7 is shown in Figure [4] and for 
rcdshift z ~ 8 in Figure [5] The data points indicate the 
maximum likelihood at each magnitude for our sam- 
ple taking into account all fields, while error bars indicate 
the smallest marginalized interval to encompass 68% of 
the likelihood for each bin. 

4.2.2. Schechter Luminosity Functions 

A determination of the ISchechterl ()1976D luminosity 
function parameters is calculated by estimating the mean 
galaxy abundance in each bin as 



1 
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t+0.5AAfu 



UV ^M uv , fc -0.5Aj\/uv 
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and then using Equations [5] and [5] to calculate the like- 
lihood of the fit parameters. The likelihoods of each 
binned sample in each field are multiplied. To im- 
prove constraints at the bright end at z ~ 7, when 



fitting a we also inclu d e dat a from the ground-based 
surveys of lOuchi et all (l2009l) . iCastellano et~aT1 pOlOl ) 
and iBowler et al.l (|2012l ). Incorporating the wide-area 
ground-based constraints is critical as even our wide area 
HST data only detects sources dimmer than Muv ~ 
21.0, or approximately 1 magnitu de brighter than M* T y 



at this redshift . As p ointed out in iRobertson I ( 2010b) 
IBouwens etaH (|2011l) . and IBradlev et all (|2012l ). there 



remains a large degeneracy between M^ v and the faint 
end slope, a without sufficient data at the bright end. 
These additional data are incorporated using the pub- 
lished data points and errors through a x 2 likelihood, as- 
suming the reported errors are Gaussian. The maximum 
likelihood Schechter function parameters are determined 
using MultiNest to conduct Bayesian inference. 

The Schechter function fit results for redshift z ~ 7 
are shown in Figure 0] and for z ~ 8 in Figure [5l 
with the maximum likelihood models shown as a black 
line. At z ~ 7, the best fit Schechter function pa- 



rameters are log 



Mi 



-3.19 



uv,*- = — 20.1 ± 0.4, and a z ^ 7 — — x. ol _ u? , 
uncertainty range for each parameter reflects the small- 
est interval in each marginalized distribution to encom- 
pass 68% of the posterior probability. At z ~ 8, the 
best fit Schechter function parameters are log 10 A 

0.35 ,„„ A ^„-3 mo „-l 
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In Figures H] and O the grey shaded 
regions denote the inner 68% of the marginalized poste- 
rior distribution in galaxy abundance at each magnitude. 

We caution the reader against an over-interpretation of 
the faintest bins in our Luminosity functions. Although 
heating of the intergalactic medium during reionization 
is expected to suppress the formation o f dwarf galaxies 
below a characteristic halo mass (e.g., iWvithe fc Loebl 
I2006t IMufioz & Loebl 1201 11 ). the density determinations 
of 4>k in our faintest bins are very sensitive both to up- 
scattering of sources below the limit and completeness 
corrections. This is largely a result of being in a regime 
where the effective volume is rapidly changing as a func- 
tion of magnitude. For example, simply dividing the 
number of observed sources, after correcting for the ex- 
pected contamination, in our faintest z ~ 7 bin by the ef- 
fective volume yields log 10 4>k ~ —2.3, much more in line 
with our best fit Schechter function parameters. Though 
we have made our best effort to quantify and account for 
corrections arising from finite size, scattering, and con- 
tamination, the situation remains difficult at the faintest 
reaches of the survey. 

4.2.3. Cosmic Variance 

Although we have not included the effects of cosmic 
variance in any of the parameter estimates given above, it 
nonetheless useful to obtain some indication of its effect. 
To first order, cosmic variance is unlikely to have a major 
effect on one of the primary goals of this paper, namely 
an estimate of the faint end slope at z ~ 7 — 8. In the 
following, we therefore give an approximate calculation 
of the effective variance arising from large scale structure. 

Density fluctuations owing to large scale modes can 
cause variations in the observed galaxy abundance be- 
yond uncertai nties arising f rom nu mber counting statis- 
tics. Following Ro bertson] (|2010bD . by using our best fit 
luminosity functions at z ~ 7 — 8, we can characterize 



these cosmic variance uncertainties for each field in our 
sample. We use the line a r pow er spectrum calculated 
with the Eisenstein & Hu ( 1998) transfer function, con- 
servatively assuming root-mean-squared density fluctua- 
tions in volumes of radius 8 /i _1 Mpc of erg = 0.9 at z = 0, 
to estimate the typical root-mean-squared density fluc- 
tuations in our survey fields at the redshifts of interest. 
To estimate the clustering bias of galaxies at these red- 
shifts, we simply match the abundance of galaxies from 
our luminosity functions with the abudance o f dark mat- 
ter halos provided by the I Tinker et al.1 ((2008) halo mass 
function, and then assign the clustering bia s of the halos 
to the galaxies assuming the bias function of I Tinker et al.l 
poll . For the UDF, to our limiting magnitude wc find 
that the typical cosmic variance (the fractional uncer- 
tainty in the total galaxy number counts owing to large 
scale structure) is acv ~ 0.30 at z ~ 7 and acv ~ 0.36 
at z ~ 8. The typical bias for galaxies in the UDF is 
b ps 5.0 at z ~ 7 and b ps 6.2 at z ~ 8. For UDF-P1, 
we find that the typical cosmic variance is acv ~ 0.33 
at z ~ 7 and acv ~ 0.38 at z ~ 8. The typical bias for 
galaxies in the UDF-P1 is b pa 5.4 at z ~ 7 and 6 ps 6.6 
at z ~ 8. For UDF-P2, we find that the typical cosmic 
variance is ctcv ~ 0.32 at z ~ 7 and crcv ~ 0.37 at 
2 ~ 8. The typical bias for galaxies in the UDF-P2 is 
b ps 5.2 at 2 - 7 and 6 ps 6.5 at z ~ 8. For ERS, we find 
that the typical cosmic variance is acv ~ 0.30 at z ~ 7 
and crcv p^ 0.34 at z ~ 8. The typical bias for galaxies 
in the ERS field is b ps 6.4 at z ~ 7 and 6 ps 7.7 at z ~ 8. 
Lastly, for CANDELS-Deep we find that the typical cos- 
mic variance is acv ~ 0.26 at z ~ 7 and trey ~ 0.30 at 
z ~ 8. The typical bias for galaxies in CANDELS-Deep 
is b ps 6.3 at z ~ 7 and 6 «s 7.6 at z ~ 8. 



5. DISCUSSION 

It is instructive to compare the Schechter function pa- 
rameters derived by our study to those of previous analy- 
ses, both at these redshifts, and below. We focus here on 
the faint end slope, a, as this is the primary parameter 
for which we expect our new UDF12 program to provide 
an improved constraint. Previous studies at redshifts 
2 < z < 6 find a remar kably consistent value of a ~ 
across this range (e.g., Oesch et alJl2007l IBouwens et al.1 
[20071 IReddv fc Steideill2b09b IMcLure et al.ll2009fl . ~ 

At the moderately larger redshift of z ~ 7, the situa- 
tion re mains much more uncertain. While lBouwens et al.l 
([20H claim a significa ntly steepened value of a = 
-2.01 ± 0.21, at z - 7, iGrazian et all (|2012f l find no 
signal of slope evolution, determining a — — 1.7 ± 0.1. 
Our determination of a z ^ — —1.871^7, though still 
consistent with —1.7, does suggest a steepening of the 
faint end slope with increasing redshift, especially when 
considering the value of a = — 1.90+ ; 15 determined by 
IMcLure et all (|2012h . who also use the new UDF12 data, 
but perform a completely independent analysis. It is re- 
assuring that the two methods agree so closely, despite 
the differing selection techniques and different potential 
sources of uncertainty. 

At z ~ 8, the existing literature largely agrees on 



a steepe ning of a, with the most recent determ ina- 
tions bv IBouwens et al.l (|2011l ). lOesch et al.1 (|2012D and 
iBradlev et al.l (|2012T ) finding values between —1.9 and 
—2.1. Extending 0.25 magnitudes fainter in UV lumi- 
nosity than any previous study, our determination of 
ot z ~8 = — 1- 94^0 24 provides increased support for this 
evolution, in concert with the a z ~s = — 2.02+g 23 found 
bvlMcLure et al l (120121) As noted by m any au thors(e.g., 
iRobertson et al.ll2010al IBouwens et al.|[2012al ). this will 
significantly increase the ability of galaxies to maintain 
the reionization of the intergalactic medium as intrinsi- 
cally faint sources become more numerous. This steepen- 
ing is also predicted by conditional luminosity function 
methods based on the evolution of the dark matter halo 
mass function (|Trenti et al.l l2010h . We also note that 
although our derived values of 4>* an d Afuv,* favor an 
decreasing <f>+ with redshift to account for the evolution 
of the luminosity function, the errors are still too large 
to rule out an evolution in the characteristic magnitude 
instead (or a combination) . We defer a significantly more 
detaile d analysis of the impact of our survey on reioniza- 
tion to IRobertson et al.l (120121). 

Along with IMcLure et all (j2Q12f > , we have uncovered 
the most comprehensive and robust sample of sublumi- 
nous high-redshift galaxies to date. At moderate magni- 
tudes, Afuv < —18.0, we achieve a more refined sample 
of dropouts, including an additional 3 (3) z-drops (Y- 
drops) not previously identified as high redshift sources 
as a result of our improved photometry. Of greater im- 
portance, though, are our advances below this UV mag- 
nitude; we discover an additional 14 sources at z ~ 7 
by virtue of our ultradeep Y105 image, as well as an ad- 
ditional 5 sources at z ~ 8, indicating the steepness of 
the faint end slope continues beyond 2 magnitudes below 
-Wuv,* a t these redshifts. Additionally, our sample is in 
excel lent agreement with the independent determination 
from lMcLure et aL (|2012t ). We note only 2 of our sources 
at z ~ 7, and 3 at z ~ 8 are not present in their final 
catalog. 

With the upcoming HST Frontier Fields observations 
scheduled to begin in Cycle 21, progress in this regime vi- 
tal to understanding if and when starforming galaxies can 
maintain reionization is sure to continue. We stress the 
gains made by UDF12 strengthen claims of an increased 
steepn ess at the faint end and, along with IMcLure et al.l 
(2012) provide a self-consistent, robust determination of 
a at redshifts 7 and 8. 
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TABLE 2 



z-drops 6.2 < z < 7.3 


ID 


RA 


Dec 


Z 850LP Yl05W Jl25W 


Jl40W 


Hieow 


References 



UDF12- 


-3746- 


-6328 


3 


:32 


:37 


.46 


-27: 


46 


:32 


.8 


28.4 ± 





.1 


27.6 


± 0.1 a 


27.5 ± 0.1 a 


27.4 ± 0.1 a 


27.3 ± 0.1 a 


1,3,9 


UDF12- 


-4256- 


-7314 


3 


:32 


:42 


.56 


-27: 


,17 


:31 


.4 


30.4 ± 





.7 


28.1 ± 0.1 a 


27, 


S 


± 0.1 a 


27.8 ± 0.1 a 


27.7 ± 0.1 a 


1,2,3,4,5,7,10 


UDF12- 


-4219- 


-6278 


3 


:32 


:42 


.19 


-27: 


46 


:27 


.8 


29.2 ± 





.2 


28.1 ± 0.1 a 


28 


.1 


± 


0.1 


28.0 


± 


0.1 


28.1 ± 0.1 a 


1,3,4,7,9,10 


UDF12- 


-3677- 


-7536 


3 


:32 


:36 


.77 


-27: 


,17 


:53 


.6 


29.0 ± 





.2 


28.2 ± 0.1 a 


28 


.2 


± 


0.1 


28.3 


± 


0.1 


28.2 ± 0.1 a 


1,2,3,4,7,9,10 


UDF12- 


-3744- 


-6513 


3 


:32 


:37 


.11 


-27: 


46 


:51 


.3 


29.5 ± 





.2 


28.4 ± 0.1 a 


28 


.3 


± 


0.1 


28.5 


± 


0.1 


28.5 


± 


0.1 


1,2,3,4,5,7,9,10 


UDF12- 


-4105- 


-7156 


3 


:32 


:41 


.05 


-27: 


47 


:15 


.6 


30.4 ± 





.7 


28.7 ± 0.1 


28 


.4 


± 


0.1 


28.3 


± 


0.1 


28.4 


± 


0.1 


1,2,3,5,10 


UDF12- 


-3958- 


-6565 


3 


:32 


:39 


.58 


-27: 


46 


:56 


.5 


29.8 ± 





.3 


28.4 ± 0.1 a 


28 


.4 


± 


0.1 


28.4 


± 


0.1 


28.4 


± 


0.1 


1,2,3,5,7,9,10 


UDF12- 


-3638- 


-7162 


3 


:32 


:36 


.38 


-27: 


47 


:16 


.2 


29.5 ± 


0. 


.2 


28.5 ± 0.1 a 


28. 


.5 


± 


0.1 


28.4 


± 


0.1 


28.5 


± 


0.1 


1,2,3,4,5,7,9,10 


UDF12- 


-4057- 


-6436 


3 


:32 


:40 


.57 


-27: 


46 


:43 


.6 


29.9 ± 





,1 


28.6 


± 


0.1 


28. 


.6 


± 


0.1 


28.5 


± 


0.1 


28.7 


± 


0.1 


1,2,3,4,5,7,9,10 


UDF12- 


-4431- 


-6452 


3 


:32 


:14 


.31 


-27: 


46 


:45 


.2 


29.6 ± 


0. 


.3 


28.7 


± 


0.1 


28. 


.6 


± 


0.1 


28.7 


± 


0.1 


28.9 


± 


0.2 


1,7,9,10 


UDF12- 


-4160- 


-7045 


3 


:32 


:41 


.60 


-27: 


47 


:04 


.5 


30.0 ± 





.5 


28.8 


± 


0.1 


28. 


.8 


± 


0.1 


28.7 


± 


0.1 


28.8 


± 


0.1 


7,9,10 


UDF12- 


-4268- 


-7073 


3 


:32 


:42 


.68 


-27: 


47 


:07 


.3 


29.9 ± 





.3 


29.1 


± 


0.1 


28. 


.9 


± 


0.1 


28.9 


± 


0.2 


28.9 


± 


0.1 


10 


UDF12- 


-3313- 


-6545 


3 


:32 


:33 


.13 


-27: 


46 


:54 


.5 


< 31 


.3 




29.3 


± 


0.1 


28. 


.9 


± 


0.1 


28.6 


± 


0.1 


29.0 


± 


0.1 


1,6,7,9,10 


UDF12- 


-3402- 


-6504 


3 


:32 


:34 


.02 


-27: 


46 


:50 


.1 


30.8 ± 


0. 


.7 


29.3 


± 


0.1 


28. 


.9 


± 


0.1 


29.0 


± 


0.1 


29.2 


± 


0.2 


7,10 


UDF12- 


-4182- 


-6112 


3 


:32 


:41 


.82 


-27: 


46 


:11 


.2 


30.3 ± 


0. 


.5 


29.3 


± 


0.1 


29. 


.1 


± 


0.1 


29.0 


± 


0.1 


28.9 


± 


0.1 


3,4,5,7,9,10 


UDF12- 


-3734- 


-7192 


3 


:32 


:37 


.34 


-27: 


47 


:19 


.2 


30.3 ± 


0. 


.5 


29.3 


± 


0.1 


29. 


.2 


± 


0.2 


29.5 


± 


0.2 


29.5 


± 


0.2 


10 


UDF12- 


-4239- 


-6243 


3 


:32 


:42 


.39 


-27: 


46 


:24 


.3 


30.0 ± 





.3 


29.2 


± 


0.1 


29. 


.3 


± 


0.2 


28.9 


± 


0.1 


29.1 


± 


0.2 


3,5,7,11 


UDF12- 


-3989- 


-6189 


3 


:32 


:39 


.89 


-27: 


46 


:18 


.9 


30.5 ± 


0. 


.8 


29.4 


± 


0.1 


29. 


.3 


± 


0.2 


29.5 


± 


0.2 


29.5 


± 


0.2 


1,10 


UDF12- 


-4068- 


-6498 


3 


:32 


:40 


.68 


-27: 


46 


:49 


.8 


30.6 ± 


0. 


.8 


29.7 


± 


0.1 


29. 


,1 


± 


0.2 


29.4 


± 


0.2 


29.6 


± 


0.2 


10 


UDF12- 


-3853- 


-7519 


3 


:32 


:38 


.53 


-27: 


47 


:51 


.9 


< 31. 


.4 




29.6 


± 


0.1 


29. 


.4 


± 


0.2 


29.4 


± 


0.2 


29.7 


± 


0.3 


2,7,10 


UDF12- 


-4472- 


-6362 


3 


:32 


:14 


.72 


-27: 


46 


:36 


.2 


30.8 ± 


1. 


.0 


29.0 


± 


0.1 


29. 


.4 


± 


0.3 


28.8 


± 


0.1 


28.9 


± 


0.1 


10 


UDF12- 


-3983- 


-6189 


3 


:32 


:39 


.83 


-27: 


46 


:18 


.9 


30.0 ± 


0. 


.1 


29.2 


± 


0.1 


29. 


.4 


± 


0.2 


29.7 


± 


0.3 


29.4 


± 


0.1 


1,11 


UDF12- 


-3736- 


-6245 


3 


:32 


:37 


.36 


-27: 


46 


:24 


.5 


30.6 ± 





.7 


29.4 


± 


0.1 


29. 


.5 


± 


0.2 


29.5 


± 


0.2 


29.6 


± 


0.2 


7,10 


UDF12- 


-3456- 


-6493 


3 


:32 


:34 


.56 


-27: 


46 


:49 


.3 


31.4 ± 


3 


.9 


29.6 


± 


0.1 


29. 


.5 


± 


0.2 


29.3 


± 


0.1 


29.2 


± 


0.1 


7,10 


UDF12- 


-3859- 


-6521 


3 


:32 


:38 


.59 


-27: 


46 


:52 


.1 


30.4 ± 


0. 


.6 


29.3 


± 


0.1 


29. 


.5 


± 


0.2 


29.9 


± 


0.3 


29.5 


± 


0.2 


7,10 


UDF12- 


-4384- 


-6311 


3 


:32 


:43 


.84 


-27: 


46 


:31 


.1 


< 30. 


.9 




29.8 


± 


0.2 


29. 


.6 


± 


0.3 


29.5 


± 


0.2 


29.9 


± 


0,1 


10 


UDF12- 


-3755- 


-6019 


3 


:32 


:37 


.55 


-27: 


46 


:01 


.9 


< 31 


.1 




29.9 


± 


0.2 


29. 


.6 


± 


0.2 


29.9 


± 


0.1 


29.6 


± 


0.2 


7,10 


UDF12- 


-3975- 


-7451 


3 


:32 


:39 


.75 


-27: 


47 


:45 


.1 


30.9 ± 


1. 


.1 


29.4 


± 


0.1 


29. 


.6 


± 


0.2 


29.2 


± 


0.2 


29.3 


± 


0.2 


7,10 


UDF12- 


-4201- 


-7074 


3 


:32 


:42 


.01 


-27: 


47 


:07 


.4 


30.6 ± 


1. 


.0 


29.6 


± 


0.1 


29. 


.6 


± 


0.2 


30.0 


± 


0,1 


29.6 


± 


0.2 


10 


UDF12- 


-4037- 


-6560 


3 


:32 


:10 


.37 


-27: 


46 


:56 


.0 


< 31. 


.3 




29.9 


± 


0.2 


29. 


.6 


± 


0.2 


30.1 


± 


0.1 


30.1 


± 


0.3 


7,10 


UDF12- 


-4426- 


-6367 


3 


:32 


:14 


.26 


-27: 


46 


:36 


.7 


30.8 ± 


1. 


.4 


29.6 


± 


0.1 


29. 


.6 


± 


0.3 


29.7 


± 


0.3 


29.5 


± 


0.2 


- 


UDF12- 


-3909- 


-6092 


3 


:32 


:39 


.09 


-27: 


46 


:09 


.2 


30.2 ± 





.6 


29.5 


± 


0.1 


29. 


.7 


± 


0.2 


30.3 


± 


0,1 


29.9 


± 


0.3 


10 


UDF12- 


-4143- 


-7041 


3 


:32 


:41 


.43 


-27: 


47 


:04 


.1 


30.6 ± 


0. 


.7 


29.9 


± 


0.2 


29. 


.8 


± 


0.2 


30.2 


± 


0.1 


30.0 


± 


0,1 


10 


UDF12- 


-3696- 


-5536 


3 


:32 


:36 


.96 


-27: 


45 


:53 


.6 


30.5 ± 





.5 


29.7 


± 


0.2 


29. 


.8 


± 


0.3 


29.3 


± 


0.2 


29.9 


± 


0.3 


7,10 


UDF12- 


-3897- 


-8116 


3 


:32 


:38 


.97 


-27: 


48 


:11 


.6 


30.7 ± 


0. 


.8 


29.9 


± 


0.2 


29. 


.8 


± 


0.3 


30.0 


± 


0.1 


29.8 


± 


0.3 


10 


UDF12- 


-4288- 


-6261 


3 


:32 


:42 


.88 


-27: 


46 


:26 


.1 


30.7 ± 


1. 





29.7 


± 


0.2 


29. 


.8 


± 


0.2 


30.0 


± 


0.3 


30.5 


± 


0.7 


10 


UDF12- 


-3817- 


-7327 


3 


:32 


:38 


.17 


-27: 


47 


:32 


.7 


< 30. 


.5 




30.0 


± 


0.2 


29. 


.8 


± 


0.3 


30.0 


± 


0,1 


30.3 


± 


0.6 


11 


UDF12- 


-4379- 


-6510 


3 


:32 


:43 


.79 


-27: 


46 


:51 


.0 


30.8 ± 


1. 


.2 


29.4 


± 


0.1 


29. 


.9 


± 


0.3 


29.8 


± 


0.3 


29.8 


± 


0.3 


10 


UDF12- 


-3691- 


-6517 


3 


:32 


:36 


.91 


-27: 


46 


:51 


.7 


< 30. 


.7 




30.0 


± 


0.2 


29. 


.9 


± 


0.3 


29.7 


± 


0.3 


30.2 


± 


0,1 


8,10 


UDF12- 


-4071- 


-7347 


3 


:32 


:40 


.71 


-27: 


47 


:34 


.7 


< 31. 


.4 




29.7 


± 


0.2 


30. 


.0 


± 


0,1 


29.4 


± 


0.2 


30.0 


± 


0.6 


7,10 


UDF12- 


-4036- 


-8022 


3 


:32 


:40 


.36 


-27: 


18 


:02 


.2 


< 31. 


.2 




30.1 


± 


0.2 


30. 


.0 


± 


0,1 


29.9 


± 


0.3 


30.7 


± 


1.3 


3,7,10 


UDF12- 


-3922- 


-6149 


3 


:32 


:39 


.22 


-27: 


46 


:14 


.9 


30.8 ± 


1. 


.0 


30.0 


± 


0.2 


30. 


.2 


± 


0.5 


30.6 


± 


0.6 


30.1 


± 


0.3 


10 


UDF12- 


-4245- 


-6534 


3 


:32 


:42 


.45 


-27: 


46 


:53 


.4 


< 30. 


.9 




30.3 


± 


0.2 


30 


.2 


± 


0,1 


29.5 


± 


0.2 


29.6 


± 


0.2 


11 


UDF12- 


-4263- 


-6416 


3 


:32 


:42 


.63 


-27: 


46 


:41 


.6 


31.1 ± 


3. 


.2 


30.0 


± 


0.2 


30 


.2 


± 


0.3 


30.3 


± 


0.5 


30.1 


± 


0,1 


10 


UDF12- 


-4090- 


-6084 


3 


:32 


:40 


.90 


-27: 


46 


:08 


.4 


< 31 


.3 




30.2 


± 


0.3 


30 


.5 


± 


0.6 


30.3 


± 


0,1 


29.8 


± 


0.3 


11 


UDF12- 


-4019- 


-6190 


3 


:32 


:40 


.19 


-27: 


46 


:19 


.0 


< 31. 


.0 




29.9 


± 


0.2 


30 


.5 


± 


0.6 


29.9 


± 


0.3 


29.8 


± 


0.2 


10 



Note. — z-drop photometry. All magnitudes listed are measured in the 70% inclusive aperture sizes listed in Section 2.4 and 
not corrected to total here. We also note that errors in magnitude space become sig nificantly asymmetric bel ow ~ IOct. We hav e 
plo tted the larger error wh ereve r appropriate. Upper limi ts are la. References arc (1) McLurc ct al. (2010) (2) Ocsch ct al. (2010), 
(3)IFinkelstein et al.l (T2010), (4) Wilkins ct al. (2011), (5) IYan et al.l 120101) (6) Lorcnzoni ct al. (2011) (7) Bo uwens et al.l (120111 )78) 
IBouwens et alj | |2011|) potential (9) IMcLure et al.l l!20Tl] ) (10) IMcLure et al.l (120121 ) robust (11) iMcLure et al.l 1 120121) potential 

a Photometric errors < 10%. 
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TABLE 3 
Y-drops 7.3 < z < 8.5 



ID 



RA 



Dec 



Yl05W 



Jl25W 



References 



UDF12- 


-3880-7072 


3 


:32 


:38 


.80 


-27: 


47 


:07. 


2 


28.0 ± 0.1 a 


27.3 ± 0.1 a 


27.3 ± 0.1 a 


27.2 ± 0.1 a 


1,2,3,4,5,7,9,10 


UDF12- 


-4470-6443 


3 


:32 


:11 


.70 


-27: 


46 


:14. 


3 


28.5 ± 


0.1 


27.8 


± 0.1 a 


27.8 ± 0.1 a 


27.8 


± 0.1 a 


1,2,3,5,7,9,10 


UDF12- 


■3952-7174 


3 


:32 


:39 


.52 


-27: 


47 


:17. 


4 


29.3 ± 


0.1 


28.3 


± 


0.1 


28.2 


± 


0.1 


28.0 ± 0.1 


1,2,3,5,7,9,10 


UDF12- 


-4314-6285 


3 


:32 


:43 


.11 


-27: 


46 


:28. 


5 


29.0 ± 


0.1 


28.3 


± 


0.1 


28.2 


± 


0.1 


28.1 ± 0.1 a 


1,2,3,4,5,7,9 


UDF12- 


-3722-8061 


3 


:32 


:37 


.22 


-27: 


48 


:06. 


1 


29.0 ± 


0.1 


28.3 


± 


0.1 


28.3 


± 


0.1 


28.3 


± 


0.1 


1,2,3,5,7,9,10 


UDF12- 


-3813-5540 


3 


:32 


:38 


.13 


-27: 


45 


:54. 





30.1 ± 


0.2 


28.6 


± 


0.1 


28.5 


± 


0.1 


28.4 


± 


0.1 


1,3,5,6,7,9,10 


UDF12- 


-3780-6001 


3 


:32 


:37 


.80 


-27: 


46 


:00. 


1 


29.8 ± 


0.2 


28.6 


± 


0.1 


28.6 


± 


0.1 


28.7 


± 


0.1 


1,3,5,6,7,10 


UDF12- 


-3764-6015 


3 


:32 


:37 


.64 


-27: 


46 


:01. 


5 


30.5 ± 


0.3 


28.9 


± 


0.1 


28.9 


± 


0.1 


28.8 


± 


0.1 


1,3,5,7,10 


UDF12- 


-3939-7040 


3 


:32 


:39 


.39 


-27: 


47 


:04. 





30.0 ± 


0.2 


29.1 


± 


0.1 


29.1 


± 


0.1 


28.9 


± 


0.1 


10 


UDF12- 


-4474-6450 


3 


:32 


:14 


.74 


-27: 


46 


:45. 





29.7 ± 


0.2 


29.0 


± 


0.1 


29.1 


± 


0.2 


29.1 


± 


0.1 


1,2,3,5,7,10 


UDF12- 


-4309-6277 


3 


:32 


:43 


.09 


-27: 


46 


:27. 


7 


30.3 ± 


0.3 


28.9 


± 


0.1 


29.4 


± 


0.2 


29.3 


± 


0.1 


1,3,5,7,11 


UDF12- 


-4309-6260 


3 


:32 


:43 


.09 


-27: 


46 


:26. 





30.1 ± 


0.2 


29.3 


± 


0.2 


28.9 


± 


0.1 


29.3 


± 


0.2 


7,10 


UDF12- 


-3463-6472 


3 


:32 


:34 


.63 


-27: 


46 


:47. 


2 


30.3 ± 


0.2 


29.6 


± 


0.2 


29.3 


± 


0.2 


29.4 


± 


0.2 


7 


UDF12- 


-3551-7443 


3 


:32 


:35 


.51 


-27: 


47 


:14. 


3 


30.7 ± 


0.4 


29.5 


± 


0.2 


29.9 


± 


0.4 


29.4 


± 


0.3 


7,11 


UDF12- 


-4336-6203 


3 


:32 


:43 


.36 


-27: 


46 


:20. 


3 


< 31. 


.8 


29.5 


± 


0.2 


29.6 


± 


0.2 


29.4 


± 


0.2 




UDF12- 


-4240-6550 


3 


:32 


:42 


.40 


-27: 


46 


:55. 





30.4 ± 


0.3 


29.4 


± 


0.1 


29.7 


± 


0.2 


29.5 


± 


0.2 


7,10 


UDF12- 


-4033-8026 


3 


:32 


:40 


.33 


-27: 


48 


:02. 


6 


30.1 ± 


0.3 


29.3 


± 


0.2 


29.6 


± 


0.3 


29.5 


± 


0.3 


3,7,10 


UDF12- 


-4308-6242 


3 


:32 


:43 


.08 


-27: 


46 


:24. 


2 


30.4 ± 


0.3 


29.7 


± 


0.3 


30.0 


± 


0.4 


29.6 


± 


0.2 


7,10 


UDF12- 


-3931-6180 


3 


:32 


:39 


.31 


-27: 


46 


:18. 





29.8 ± 


0.2 


29.3 


± 


0.2 


29.3 


± 


0.1 


29.6 


± 


0.2 


10 


UDF12- 


-3934-7256 


3 


:32 


:39 


.34 


-27: 


47 


:25. 


6 


31.2 ± 


0.7 


29.8 


± 


0.3 


30.0 


± 


0.3 


29.7 


± 


0.2 


7,11 


UDF12- 


-3881-6343 


3 


:32 


:38 


.81 


-27: 


46 


:34. 


3 


30.6 ± 


0.1 


29.9 


± 


0.1 


30.0 


± 


0.4 


29.8 


± 


0.3 


11 


UDF12- 


-3681-6421 


3 


:32 


:36 


.81 


-27: 


46 


:42. 


1 


< 31. 


.4 


30.0 


± 


0.1 


30.0 


± 


0.4 


29.8 


± 


0.3 


7,11 


UDF12- 


-4294-6560 


3 


:32 


:42 


.94 


-27: 


46 


:56. 





30.8 ± 


0.6 


29.9 


± 


0.3 


30.5 


± 


0.7 


29.9 


± 


0.3 


11 


UDF12- 


-3920-6322 


3 


:32 


:39 


.20 


-27: 


46 


:32. 


2 


< 31. 


.9 


29.9 


± 


0.3 


29.7 


± 


0.3 


29.9 


± 


0.3 


10 


UDF12- 


-3858-6150 


3 


:32 


:38 


.58 


-27: 


46 


:15. 





30.3 ± 


0.3 


29.8 


± 


0.3 


30.3 


± 


0.4 


29.9 


± 


0.3 


11 


UDF12- 


-4344-6547 


3 


:32 


:43 


.44 


-27: 


46 


:54. 


7 


< 31. 


.9 


30.0 


± 


0.3 


30.1 


± 


0.4 


30.0 


± 


0.3 


10 


UDF12- 


-4101-7216 


3 


:32 


:41 


.01 


-27: 


47 


:21. 


6 


30.4 ± 


0.2 


29.9 


± 


0.3 


30.0 


± 


0.3 


30.9 


± 


0.8 


10 



Note. — Y-drop photometry. All magnitudes listed are measured in the 70% inclusive aperture sizes listed in Section 
2.4 and not corrected to total here. We also note that errors in magnitude space become significantl y asymmetric below ~ 
10< 7. We have plotted the l arger error whereve r appropriat e. Upper limi t s are la. References a re ( DIMcLure et al.1 fWuj) 
(2llOesch etaLl 120101), (3) IFinkelstein et all 120101) (4llWilkins et al.l 1 12011), (5)IYanet al.l 112010) (6) ILorenzoni eTaU 
(2011) (7) Bou wens et all 1120111) (8) IBouwens et al.1 1(20111 ) potential O MMcLure et al.1 l |201ll ) (lO HMcLure et al.l l|2012l ) 
robust (ll) IMcLure et al.l 1120121 ) potential 

a Photometric errors < 10%. 



TABLE 4 
SWML determination of the 

z ~ 7 LF 



Af uv log 4> k [Mpc" 3 



mag x \ 



-20.65 
-20.15 
-19.65 
-19.15 
-18.65 
-18.15 
-17.65 
-17.15 



4.29 
3.71 
3.31 
3.02 
2.98 
2.56 
2.23 
3.03 



.+0.29 




'-2.34 
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TABLE 5 
SWML determination of the 

z ~ 8 LF 

M uv log 4> k [Mpc- 3 mag" 1 ] 

-22.00 < 5.01 

-21.50 -5.02±°;^ 

-21.00 -4.28I024 

-20.50 -4.15i{j;43 

-20.00 -3- 54 ±o.06 

-19.50 -3.34tp|f 

-19.00 - 2 - 97 -o.20 

-18.50 -2.9ll°;24 

-18.00 -2.6ll[5:2o 

-17.50 -2.57ip:?4 



5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 



5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 





0.2 



Jl25 ■ Hi. 




Fig. 1. — Left: Number of 2850-drop contaminants per UDF12 field as a function of the color cuts in 2350 — Vios and F105 — Ji25- The 
selection criteria are defined such that the 2850 — F105 must be greater than the value on the y-axis to be selected, and theYios — J125 less 
than the value on the x-axis. Right: As left, but for Yios-drops. Our chosen cuts are marked by the black 'x' in each figure. We refer the 
reader to Figure \3\ for a visualization of our final cuts. 



The UV Luminosity Function of Star-Forming Galaxies at z ~ 7 and 8 
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Fig. 2. — Left: Selection function for zsso-drop galaxies in the UDF, as a function of Mtjv and redshift, constructed using the simulations 
discussed in section 4.1. Even at bright magnitudes, maximum efficiency is only ~0.65 due to the area subtended by other objects and 
our strict optical non-detection criteria, which result in a small fraction of true high-redshift sources being excluded. Right: Equivalent 
selection function for Yios-drops. 



UDF12 z-drops 




UDF12 Y-drops 



(3 = -2.0 
P = -3.0 




Jl25 Hj, 



Fig. 3. — Color-color diagram of galaxies selected as 2850-drops (left) and Yios-drops (right). Red and orange curves show the tracks of 
synthetic high-redshift galaxies for various UV continuum slopes /3. Light blue points denote 1-cr upper limits. 
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Fig. 4. — The luminosity function of star-forming galaxies at z ~ 7 from the zsso-drop sample. Black points were determined using 
the UDF12 data set and other HST data mentioned in this work. Red points denote wide area ground based data increase the range in 
luminosity. The black line defines the maximum likelihood Schechter luminosit y function and the shaded grey region denotes the 68% 
confidence interval. The green dashed line denotes the fit of|McLurc ct al. ( 2013). 



The UV Luminosity Function of Star-Forming Galaxies at z ~ 7 and 8 
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Schenker et al. <p{m) 

$ Schenker et al. SWML 

- - McLure et al. (2013) 
f Bouwens et al. (2011) 
f Bradley et al. (2012) 



-23 -22 -21 -20 -19 -18 
m [AB Magnitude] 



-17 



-16 



Fig. 5. — The luminosity function of star-forming galaxies at z ~ 8 from the lios-drop samp le. Black points were determined using the 
UDF12 data set and other HST data mentioned in this work. Red points denote data from the Bradley et al. (2012) analysis of the BoRG 
fields that increase the range in luminosity. The black line defines the maximum likelihood Schechter luminosity function and the shaded 
grey region denotes the 68% confidence interval. The green line denotes our fit when removing the BoRG datapont s; we note that our fi t 
to the faint end slope is remarkably insensitive to their inclusion/exclusion. The green dashed line denotes the fit of McLure et al.l H2012D . 



